%%%% Clean all results from SE run
clear

load resultsSE.mat

trials = length(hv);

% Reports Covariance results

for i=1:trials
    sse      = sum(seresults(i).grad.^2);
    varcovar = gprime * inv(seresults(i).hess2) * gprime' / obs;
    sevec    = sqrt(diag(varcovar));
    
    col(:,i) = [seresults(i).hv;seresults(i).grad';sse;sevec];
end

col = fliplr(col);